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Abstract 

In this paper we introduce an optical approximation into the theory of impedance calculation, one 
valid in the limit of high frequencies. This approximation neglects diffraction effects in the radiation 
process, and is conceptually equivalent to the approximation of geometric optics in electromagnetic 
theory. Using this approximation, we derive equations for the longitudinal impedance for arbitrary 
offsets, with respect to a reference orbit, of source and test particles. With the help of the Panofsky- 
Wenzel theorem we also obtain expressions for the transverse impedance (also for arbitrary offsets) . 
We further simplify these expressions for the case of the small offsets that are typical for practical 
applications. Our final expressions for the impedance, in the general case, involve two dimensional 
integrals over various cross-sections of the transition. We further demonstrate, for several known 
axisymmetric examples, how our method is applied to the calculation of impedances. Finally, we 
discuss the accuracy of the optical approximation and its relation to the diffraction regime in the 
theory of impedance. 
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I. INTRODUCTION 



The calculation of the impedance for the elements of a vacuum chamber system and 
the associated calculation of beam dynamics effects, such as beam instabilities or wakefield 
induced omittance growth, are important elements in the design of a modern accelerator. 
Sophisticated computer programs are routinely used for such calculations, and in many cases 
they can successfully treat complicated geometries, like those found in real vacuum systems. 
There are, however, cases where the simulations are pushed to their limits in resource 
requirements of both memory and cpu time, and in their ability to yield an accurate result. 

One example is the case of very short bunches, like those envisioned in future linear 
colliders and future light sources. For example, the final bunch length in the current design of 
the International Linear Collider (ILC) [T] is 300 microns; the final (rms) bunch length in the 
Linac Coherent Light Source [5j is 20 microns and in the European XFEL [3] it is 25 microns. 
Since numerical calculation of the short-range wake requires a spatial mesh size equal to a 
fraction of the bunch length, submillimetcr bunches represent a challenging computational 
task. Another example where direct numerical calculation is difficult is related to long, small- 
angle tapers which are often used to minimize the abruptness of vacuum chamber transitions. 
For example, collimators with such tapers will be used in the post-linac collimation section of 
the ILC. In such cases numerical solution of Maxwell's equations requires a large number of 
mesh points to fully cover the length of the transition or collimator. The difficulty becomes 
especially pronounced for short bunches. 

The difficulty in both examples mentioned above is associated with a small parameter. 
For short bunches such a small parameter is the ratio of bunch length a z to typical size b 
of the structure (in the vacuum chamber) that generates the impedance. If we denote by 
A the inverse wavenumber c/lu (c is the speed of light and to the characteristic frequency 
of interest) then A ~ cr z , and the small parameter for the problem is A/6. It has long 
been known that effective utilization of this small parameter may allow one to simplify 
the impedance problem, and several analytical results are available in the literature for the 
impedance at high-frequencies. They include the impedance of a step transition |U [5] and 
the diffraction model for the impedance of a cylindrical pillbox cavity [6j |7j [H [9]. More 
recently, a parabolic equation method was developed that provides a simplified treatment of 
diffraction effects at high-frequencies [10] . 

In the case of diffraction theory, the calculation takes into account the fact that radiated 
electromagnetic fields do not propagate along straight lines. A Fresnel type integral from 
the diffraction theory of light is used to evaluate the electromagnetic energy that enters 
into the cavity region [TJ. This energy is associated with the energy lost by the beam 
and is thus related to the real part of the impedance. In this paper we will show that in 
many cases the same kind of argument can be applied to the calculation of impedance in 
an approximation that we call the optical approximation (or optical regime) in the theory 
of impedance. In this approximation we assume that the electromagnetic fields carried by 
a short bunch propagate along straight lines equivalent to rays in the geometric optics. 
An obstacle inside the beam pipe can intercept the rays and reflect them away from their 
original direction. The energy in the reflected rays is associated with the energy radiated 
by the beam, which can then be related to the impedance. Note that this kind of argument 
has been used in the past in the case of step-in and step-out transitions in a round pipe, 
where the impedance was related to the energy "clipped away" from the beam by the step 
[51 [TT]. In a recent paper this approach has been extended and applied to the calculation of 
impedance for 3D collimators [T2"] . 

If I is the length of an obstacle and b is the minimal cross-section size of the beam pipe, 
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the conditions for the optical approximation are 

b 2 

A<&, (1) 

A 

The first of these two conditions requires the size of the obstacle be much larger than the 
reduced wavelength of the radiation. The right hand side of the second inequality has a 
meaning of the length over which diffraction effects become significant, and this relation 
guarantees that such effects give only a small correction to those of geometric optics. Note 
that even a small-angle taper of angle 8 ~ b/l can be described in the optical approximation 
for short enough bunches, if X/b <C 9. The quantity b 2 /X can also be interpreted as a catch- 
up distance over which radiation, generated by the head of a beam, reflects from a side wall 
of radius b and reaches the beam tail at length a z ~ A behind the head. Thus the second 
condition of Eq. for the applicability of the optical approximation is that the object is 
short compared to the catch-up distance. 

An analogous problem in geometric optics would be a body with transverse size b and 
length I illuminated by light, see Fig. [l] The energy reflected by the body can be simply 

IKE 

/ 

FIG. 1: In geometric optics, light incident on a metallic object from the left is reflected by the 
surface creating a shadow behind the object. The reflected light can be considered as radiation 
emitted by the object due to currents induced in the metal by the incident electromagnetic field. 

calculated, in the optical approximation, as the energy incident upon the cross-sectional 
area of the body. Such calculation, however, is only valid if the length of the obstacle I and 
the transverse size b obey the second condition in Eq. In the opposite limit (I 3> b 2 /X), 
diffraction effects become dominant. 

In this paper we assume perfect conductivity in the vacuum chamber wall. Our goal is 
to justify the optical approximation at high-frequencies and to demonstrate how one can 
apply it to calculations of both the longitudinal and transverse geometric impedances. We 
take two approaches to the problem. The first one is based on the general energy balance 
equation in electromagnetic theory, which relates the longitudinal impedance to the energy 
radiated from a transition in a beam pipe. In the past, calculations of impedance based 
on this relation were carried out by one of the authors in axisymmetric |13j as well as 
rectangular |14j geometries. The derivation is similar to the approach developed in Ref. |15j 
for the problem of laser acceleration in vacuum. As it turns out, in the most general case 
of unequal offsets of the leading and trailing particles, this approach gives the sum of the 
longitudinal impedances symmetrized over the coordinates of the particles, which does not 
provide enough information to obtain the transverse impedance for the transition. Our 
second approach uses a so-called indirect integration method developed in Refs. [TSJ 02] ■ 
Although not as physically transparent as the energy balance method, this method — with 
some reasonable assumptions regarding the formation length of the wake field — gives a 
simple expression for the longitudinal impedance in the most general case. 
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We want to emphasize here that although our result is not derived formally from first 
principles, it is based on a combination of exact consequences of Maxwell's equations and 
simple physical arguments that follow from the geometric optics. Our result is applicable to 
an arbitrary three dimensional transition and allows the incoming and outgoing beam pipes 
to have different cross-sections. Practically important examples of impedance in the optical 
regime are considered in a companion paper |18j where we also make a detailed comparison 
with computer simulations and find excellent agreement between the theory and simulations. 

This report is organized as follows. In Section[TT]we derive equations that relate the longi- 
tudinal impedance to the integrals of the Fourier transformed Poynting vector over surfaces 
located far from the transition. In Section [nT| we calculate the contribution to the impedance 
of the static fields carried by the charges while they are in the incoming and outgoing pipes 
and far from the transition. The derivation in these two sections is general and does not 
make any high-frequency assumptions. In Section IV the contribution to the impedance of 
the radiation field is computed in the optical approximation. In Section [V] we invoke the 
indirect integration method to find the wake field for unequal offsets of the particles. In 
Section |VI[ using the Panofsky-Wenzel theorem, we derive expressions for the transverse 
impedance in the optical approximation. Although our results for the transverse impedance 
are applicable for arbitrarily large transverse offsets, we further specialize them for the case 
of the small offsets usually assumed in practical applications. Example impedance calcu- 
lations are given in Section VII where we compute both the longitudinal and transverse 
impedances of a short, round collimator, and round step-in and step-out transitions. We 
show that our results for these cases agree with known expressions found in the literature. In 
Section [VlII| we discuss the relation between the optical regime and diffraction theory, and 
establish the accuracy of the optical approximation. Finally, in Section IX we summarize 
the main results of the paper. 

We use Gaussian system of units throughout this paper; to convert impedances to MKS 
units, one multiplies them by the factor Zqc/Att, with Zq = 377 £1. 



II. LONGITUDINAL IMPEDANCE AND THE ENERGY RADIATED BY 
CHARGED PARTICLES 

It is well known that the real part of the longitudinal impedance is related to the energy 
loss of the beam [T5]. Since we assume perfect conductivity in the metallic walls, the losses 
are due to the beam-induced radiation that propagates away from the transition. In this 
section we derive a formula that relates the impedance to integrals of the field over remote 
surfaces far from the transition, and which is valid for a transition of general shape and 
incoming and outgoing beam pipes of arbitrary cross-section. 

We consider transitions between two beam pipes of cross-sections Sa and Sb, respectively, 
as shown in Fig. [2] Sa and Sb are arbitrary, other than their axes are parallel to each other 
and also parallel to the so-called design orbit of the beam. We let the design orbit, in turn, 
define the z-axis. We assume that a charge qi moves at offset V\ = (a: 1,2/1) with respect to 
the design orbit, and a charge (72 moves at offset r 2 = (%2, 2/2) and at distance s behind (for 
s > 0) the leading charge; both charges are assumed to move at the speed of light c. We 
want to calculate the longitudinal wake Wii(ri,r2, s) given by 

c f°° 

w»(r 1 ,r 2 ,s) = / Ei jZ (r 2 ,z= ct - s,t)dt, (2) 

where E\ tZ (r, z, t) is the z-component of the electric field generated by the leading particle 
at the position of the trailing one. The longitudinal impedance is related to the wake by 
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FIG. 2: Geometry of a transition from pipe A to pipe B. The source and test particles move from 
left to right along the z-axis. The aperture area 5 ap is the minimal cross-section that connects the 
two pipes. 

the Fourier transform 

1 f°° 

Z||(ri,r 2) w) = - / dsw ll (r u r 2 ,s)e ius / c . (3) 
c Jo 

Let us denote the electric and magnetic fields of the leading charge by E\[r,z, t) and 
H\(r, z,t), and the electric and magnetic fields of the trailing charge by E 2 (r, z, t) and 
H 2 (r,z,t). The field (-E7i,i3i) is due to the charge density of the first particle, pi — 
q\8{z — ct)S(r — ri), and the field (E 2 ,H 2 ) is generated by the charge density of the second 
particle, p 2 = q 2 S(z—s—ct)S(r—r 2 ). The total field is the sum E = Ei+E 2 , H = Hi + H 2 . 

Using the energy balance equation in electrodynamics [20] we can calculate the energy 
lost by both charges by integrating the Poynting vector over remote boundaries Sa and Sb 
located, respectively, far to the left of the incoming pipe and far to the right of the outgoing 
pipe 



W = ~ / dt [ n-(Ex H)dS . (4) 

' '-co JSa+Sb 
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Here n is the unit vector perpendicular to the surface area, and we use the short-hand 
notation J Sa+Sb for the sum of the integrals over Sa and Sb- In this equation the vector n 
is oriented in the outward direction — it is parallel (antiparallel) to the z axis on Sb (Sa)- A 
positive value of W means that the charges lose energy. The contribution of such integrals 
over the metallic surface of the pipes and the transition vanishes because we assume that 
there are no losses in the wall. It follows from the energy balance equation that the radiated 
energy given by Eq. Q is equal to minus the energy change of both particles 

/CO /"CO 
E z (n,z = ct,t)dt — cq 2 / E z (r 2 , z = ct — s,t)dt . (5) 
-co J —oo 

The first term in this equation is the energy loss of the leading particle, and the second 
term is the energy loss for the trailing one. Note that strictly speaking the integration in 
Eq. ([5]) should only be performed over the region between the surfaces Sa and Sb] however, 
assuming that the surfaces are located sufficiently far from the transition (to where the 
interaction of the charges vanishes) , we can extend the limits of integration to infinity. 

The electric field at the location of the leading particle is equal to E\, z {r\,z — ct,t) + 
E 2>z {ri, z = ct, t), and the electric field at the location of the trailing particle is Ei :Z (r 2 , z — 
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ct — s,t) + E 2jZ (r 2 , z = ct — s,t). We can cast the formula for W as follows: 

W = qfUt + qi q 2 U 2 + q 2 2 U 3 , (6) 



where 



c f°° 

Ui = — - / E hz {r u z = ct,t)dt 

qi J-oo 

-f 

3i J-< 

/CO 
E 2tZ {r 2 ,z = ct-s,t)dt. (7) 
-oo 



U 2 = -—j E ltZ (r 2 ,z = ct - s,t)dt - — I E 2 ^ z (r 2 , z = ct,t)dt 



c 
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Since we assume that the particles move at the speed of light, due to causality the trailing 
particle does not affect the leading one. This means that E\^ z (r 2 , z = ct — s,t) is not equal 
to zero only if s > 0, and, similarly, E 2 , z (r 2 ,z = ct,t) does not vanish only for s < 0. This 
observation allows us to rewrite U 2 as 

U 2 = h(s)wi\(n,r 2 ,s) + h(-s)wi\(r 2 ,n,-s) , (8) 

where h(s) is the unit step function. 

The next step is to take the Fourier transform in time, introducing the fields £ and Ti., 



£(r,z,w) 
7i(r, z, lo) 



J 2nJ_ oc \ H(r,z,t) j 



(9) 



Similarly, we also define the Fourier components £i(r, z, u>) and Hi(r, z,u>) for the elec- 
tromagnetic field of the first charge. Because the fields E and H are real, we have the 
relations 

£(r,z,-u) = £*(r,z,u) , H(r,z,-u) =H*{r,z,u) , (10) 

with equivalent relations for £\ and Ti,\. The asterisk in these equations denote complex 
conjugation. As for the Fourier components of the second charge, we define them in a way 
that explicitly separates the phase factor introduced by the distance s between the particles: 

£ 2 (W) 1 = 1 Wc r ^ f E 2 (r,z,t) X 
H 2 (r,z,(j) J 2tt 7_ oc \ H 2 {r,z,t) J v ; 

We then have £ = £ x + e^ s l c £ 2 and H = Hi + e luJS / c H 2 . 

From the linearity of Maxwell's equations it follows that the electromagnetic field (£ 1 ,'H.\) 
is generated by the charge density pi equal to the Fourier image of p\ 

Pi = ^J dt Pl e^ = ^-e^6(r - n) , (12) 

and the field (£ 2 ,H 2 ) is generated by the charge density 

1 f°° a 

h = —e-^ c / dtp 2 e iut = ^ e ^ z ' c 5(r - r 2 ) . (13) 
27T ./_„ lire 
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Note that because of the extra phase factor e~ ius ' c in the definition (111 the Fourier com- 
ponents p± and p 2 are now "in phase". Since p 2 has the property p 2 (io) — p 2 {—to), the fields 
(£ 2 ,7i 2 ) satisfy the same relations as Eqs. (10 1. 

Using Parseval's theorem for Fourier transforms we can express W in Eq. ^ in terms of 
£{lo) and H(u): 



W = 



du) 



[£{w) x H*(u>)\ ■ ndS, 



(14) 



'Sa+Sb 

or, equivalcntly, in terms of the fields of the first and second particles 



W 



du> 



du> 



du> 



Sa+Sb 



Sa+Se 



Sa+Sb 



[fi(w) x • ndS 

" s/c £i{w) x H.* 2 {ia) 
[£ 2 (uj) x n* 2 (uj)} ■ ndS. 



e luJs/c £ 2 (to) x ni(u) 



■ ndS 



(15) 



We will now prove that each of the three terms on the right hand side of Eq. (|6| is equal 
to the corresponding term in Eq. (151. The proof is based on the linearity of Maxwell's 
equations and the fact that expressions ^ and ( 15 1 are equal for arbitrary values of q\ and 
Q2- If 1i — 0, then there are only the first terms in these two equations, hence they are 
equal. If 171 = 0, then there are only the third terms, and they are also equal. Hence the 
second terms must be equal, too. Using Eq. (|8| we obtain 



h(s)w\\ (ri, r 2 , s) + h(-s)w\\ (r 2 , n, -s) = 



2qi<72 

c 

2gi<72 J- 



du> 



duj i 



Sa+Sb 
iujs j c 



Sa+Sb 



[£i(lu) x H 2 (w) - 



e ibJs/c £ 2 (u) x nl(uj) ■ ndS 
£* 2 (oj) x Hi(w)] ■ ndS. 



(16) 



In the last integral we changed the integration variable lu — ► —uj in the second term of 
the integrand and used the relations £ 2 {— uj) = £ 2 {oj) and ?"£*(— ui) = T~t\{u)). Taking the 
Fourier transform of this equation and using Eq. ([3]) gives 



Z(r\,r 2 ,oS) = Z\\(n,r2,w) + Z\\(r2,ri,w 

7TC 

Q1Q2 JS A +S E 



[£\{ijS) x 7*2(0;) + £* 2 {uj) x Wi(oj)] ■ ndS . 



(17) 



The quantity Z in this equation is the sum of the impedances symmetrized with respect to 
the offsets of the leading and trailing particles. If r\ — r 2 , then the left hand side is equal 
to 2Zii(ri,ri,u>) and this equation gives the longitudinal impedance for particles that have 
the same offset. However, for the general case of unequal offsets, the equation gives only 
the sum of the impedances Z. 

Note that up to this point we did not use any approximation in our derivation — our result 
is valid not only in the optical regime, but for arbitrary frequency lu. Eq. (17), relates the 
sum of two longitudinal impedances to the interference term between the electromagnetic 
fields of charges 1 and 2 in the energy flow through the remote boundaries Sa and Sb- In 
the next sections we will split the right hand side of Eq. (17 1 into several distinct parts and 
then calculate the parts individually in the optical regime. 
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III. CONTRIBUTION TO IMPEDANCE OF STATIC AND RADIATION FIELDS 



There are several contributions to the integral in Eq. (17). The first one comes from 



the field that is carried by particles when they cross the surface Sa on their way to the 
transition, and the second one arises when they pass through the surface Sb after leaving 
the transition. We call these fields the static fields because they do not vary in time in the 
beam frame and they can be calculated from a time independent system of equations after 
a trivial change of variables z — ct — > £ in the laboratory frame. We emphasize here that 
these are not radiation fields; nevertheless, they need to be included into the original energy 
balance in Eq. Q. 

We begin calculation of the static field contributions from pipe A. The electric fields £i,a 
and £ 2 , A i n the straight pipe can be represented in terms of the potential function <fi, 



£\.A = 



9l juz/c 
27TC 



£2, A = 



9j ^J.ujz/c 

~2ttc 



V02,A(r) , 



(18) 



where V = x d/dx + y d/dy. The potential <fi satisfies Poisson's equation with source terms 
given by the charge densities of Eqs. ( 12 ) and ( 13 1: 



V 2 0i !j4 (r) = -4irS(r - n) , V 2 ^ 2 ,a(t") = -4vr£(r - r 2 ) 



(19) 



with boundary conditions <j>\ a — 4*2, A = on the wall of pipe A. The magnetic field in pipe 
A is: Hi a = z x £\ a and H2 a = z x £ 2 a- We denote by Za the contribution of these 



fields, when integrated over the surface Sa, to the integral Eq. (17) 
2ttc 



9i92 Js A 



f £i(uj) ■ £ 2 (uj)dS = f X7(f) 1A -X7(f) 2 A dS. 

Js A ~ 2ttc J Sa 



(20) 



The minus sign in this equation is due to the fact that the normal vector to Sa is oriented 
in the direction opposite to the z axis. 



In the same way we calculate the contribution Zb to Eq. ( 17 1 from the remote boundary 
Sb- The electric field on this surface £\ t B and £ 2 ^b is given by the same equations ( [18} and 
(19 1 with the index A replaced by B, and the boundary conditions for 4>i,b and 02,fl being 



zero on the metallic surface of beam pipe B. However, because the direction of the normal 
vector to Sb is along the z axis, Zb has a positive sign 



1 

2nc 



V01 B ■ V0 2 ,B dS . 



(21) 



Note that, strictly speaking, the integrals in Eqs. (20) and (21) diverge because the field 
has a singularity at the location of charges q\ and q 2 . However in the sum Z A + Zb the 
singular contributions cancel, and the result is finite. 

In addition to the terms Za and Zb there will be a contribution to the impedance due to 
the radiation field emitted in the transition region. We will denote this contribution i? ra d, 
so that 



Z = Z A + Z B + Z rad . (22) 
The quantity Z ra( j will be evaluated in the next section using the optical approximation. 
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IV. RADIATION FROM THE TRANSITION IN THE OPTICAL APPROXIMA- 
TION 



In the optical regime the radiated energy has several terms. The first term is the energy 
that is radiated by reflection from the narrowing part of the pipe in the aperture area S a p 
(see Fig. |2j. We denote the cross-section of the transition that complements S ap to Sa by 
Sa — S ap , and, similarly, the cross-section that complements S ap to Sb by Sb — S ap . The 
incident energy flux within this area will be "clipped away" from the beam and converted 
into a radiation field. This radiation may go in the backward direction if the narrowing 
region is steep enough (e.g., a diaphragm, or a 90-degree, abrupt step-in pipe radius), or, 
for a small-angle taper, it may go in the forward direction. Since the incident energy is the 
static field in pipe A, the contribution Z la ^ i of the clipped energy is given by the same 



Eq. (20), but with the integration now carried out over the "clipping" area Sa — S ap , 

Z Iad ,l = 7T~ I Wi.A ■ dS . (23) 

The sign in this equation is positive because the radiation propagates from the transition 
region to infinity. 

Through the aperture S ap connecting pipes A and B, charges carrying the static potential 
4>a enter into pipe B (for brevity, we momentarily drop indices I and 2). The fields will 
eventually change in such a way that at a large distance from the transition the particles 
will carry the potential 4>b- In the course of this restructuring, there will be additional 
radiation emitted in the forward direction. We denote the contribution of this radiation to 
the impedance by Z ra d,2 and calculate it in the following way. We represent the potential 
of the charges immediately after passing through the aperture as a sum of the potential 4>b 
occupying area Sb, the potential — 4>b over the area Sb — S ap , and the potential 4>a — <t>B 
over the area S'ap. From the linearity of Maxwell's equations, the first field proceeds with 
the charges as a new static field in pipe B, and the last two transform into radiation. The 



energy integral of the interference fields carrying this radiation is similar to Eqs. (20) and 
(21 1 with the difference that over the area S'ap one needs to use (f>A — 4>b, and outside of it, 



the potential ~(\>b 



Z la d,2 = J— [ V(4>1,A ~ 01, B) ■ V(02,A - 4>2,B) dS + — *— / V(j>l,B ■ V02.S dS . 

27TC Js ap 2ttc Js B -s ap 

(24) 

A slightly different view of the derivation of Eqs. (23 1 and (24 1 is illustrated by Fig. [3] 
According to this view, since the length of the transition is shortfto satisfy Eqs. (jlj], it can 
be treated as an infinitely short protrusion located at z = 0. The incident electromagnetic 
field in pipe A, given by potential cf>A, generates a radiation field which will be reflected back 
into pipe A. The radiation field satisfies the boundary conditions <j> — —(f) a on the left side of 
the protrusion, at z = —0, indicated in Fig. [3] This condition follows from the requirement 
that the tangential electric field on the surface of the metal is zero. The boundary conditions 
on the side of pipe B, at z = +0, are also indicated in the figure. The radiation field that 



1 For a 3D transition with a complicated geometry the general rule for finding the aperture area S'ap is the 
following. Assume that the transition is illuminated by a set of parallel rays of light that propagate from 
pipe A along the z axis. Then the cross-section of the illuminated area in pipe B gives S'ap- 
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A <t> = o 
<l>=-<i> A 



* = Mb B 



(j> = 4>j 



■ ,. 

z 



FIG. 3: The transition of Fig. [2] represented schematically as an infinitely thin protrusion. Also 
given are the boundary conditions for the potentials connected to the radiation fields on both the 
A and B sides of the protrusion. 



satisfies these boundary conditions, when summed with the static field corresponding to 
potential <pg, gives a field equal to zero on the metallic surface of the protrusion, and equal 
to 4>A over the aperture S ap . The interference terms between the radiation fields of the first 



and the second charges lead to the integrals Eq. ( 23 1 and (24 1 on the left and the right side 
of the transitions, respectively. 

Collecting now all contributions in Eqs. (20 1, (21 1, (23 1 and (24 1 we arrive at the final 
result 

Z = — f V&,b ■ V0 2 , s dS- — I (V0i. A • V0 2 , B + V0i, B • V0 2 . A ) dS , (25) 

JSb 2lTC Js ap 

where the first integral is carried out over the cross-section of pipe B, and the second integral 
over aperture >S a p. Note that the impedance given by Eq. (25) is real and does not depend 
on frequency to. 

As was emphasized above, the energy balance equation which was used to derive Eq. 
(25 1 gives the summed impedance Zn (n, r%, uS) = Zy(r2,ri,w) + Zn (n, T2, u>), which is 
symmetrized over the offsets of the leading and trailing particles. This suffices to give us 
the longitudinal impedance for the case when leading and trailing particles follow the same 
path. However, for the transverse impedance (which we will discuss below) we will need 
to first find Zn (r\, r 2 , ui) alone (Z\\ is not sufficient). To find Zn (ri, r^, oS) in the optical 
approximation, we will, in the next section, invoke a more formal approach based on a 
so-called indirect integration method. 



V. INDIRECT INTEGRATION METHOD 



The indirect integration method developed in Refs. [TH1 HZ] reduces the calculation of the 
wake field to the integration along a straight path up to some point zq located in the exit 
pipe B. The contribution from the remaining path z > zq is expressed through the solution 
of an auxiliary problem at the cross section z — z . 

By choosing zq in pipe B immediately after the transition, we note that, in the optical 
regime, the contribution of path z < zq can be neglected. Indeed, the wake accumulates 
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over the catch up distance ~ b 2 /X, which, according to Eqs. ([!]), is much larger than the 
transition length I. This observation greatly simplifies the problem and allows us to compute 
the impedance from the part of path with z > Zq only. This contribution is |16j 



w(n,r 2 ,s) 



where $ satisfies the equation 
V 2 $(r,s) = 



d_ 

dz 



at 



Er z (r,z = 0,t = s/c) , 



(26) 



(27) 



with the boundary condition $ = on the wall of pipe B. The quantity El° is (according 
to the terminology of Ref. [16]) the "scattered" electric field of the leading charge — it is 
obtained by subtracting from the total field of this charge its static field in pipe B. In Eq. 
(27 1 we set zq = and suppress the argument r% in the function <I> — the dependence on this 
argument is clear from Eq. (27 1 where the electric field Ei iZ is generated by the leading 
particle moving with offset r\. 

We now transform Eq. (27) in a way that eliminates the longitudinal component E z . 

^dEf/dt = (V x H sc ) z we obtain 



Using dEf/dz + V 



E s l = and c 



dE\- z 
dz 

at 



= -V • E 



(V x Hf ± ) z 



(28) 



where the symbol _L refers to the components of the field perpendicular to the z axis (we 
also remind the reader that V is a two-dimensional operator in the x — y plane) . 

Because we now use the time domain representation for the fields, we need to take the 
inverse Fourier transform of Eq. (T8| (and of a similar equation for pipe B) to find the static 
fields of the particle. This gives for E\ 



E lfA (r,z,t) 



6(t-z/c)V<t> lt A(r). 



E 1B (r,z,t) 



6(t-z/c)V<Px, B (r), (29) 



with the magnetic fields given by Hi a — * x El,A an d Ht g = z x E lyB . 

The crucial step in the derivation is to notice that, in the optical regime, after passage 
through the transition region, the static field of particle 1 is "scraped off" outside the 
aperture S'ap. The field left with the charge is equal to E\ t A but only within the area of 
S'ap (the field in S B — S ap is zero) . To find the scattered field we need to subtract from the 
"truncated" field E\ a the static field E% B of charge 1 in pipe B, which gives 



El c ± (r,z,t) = 



-(?i/c)Vfo M (r) - cj> hB (r)}5(t - z/c) 



-{qi/c)V<h., B {r)6{t - z/c) 
with the corresponding magnetic field given by 



m o B 



SR — S; 



ap j 



H 



1_L 



z x Ei 



(30) 



(31) 



We substitute Eqs. (30) and (31) into Eq. (28 1 and use the result as the right hand side in 
Eq. ([27} for It follows from Eq. pi} that(V x H\ c ± ) z = V • El c ± , which gives 



V 2 $(r, s) = -2V • Ef ± (r, z = 0,t= s/c) . 



(32) 
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This equation is solved by first noting that — (47t) _1 (/>b(t") is the Green function in Sb (4>b 
satisfies Eqs. (19 1 in which A is replaced by B) and then using the symmetry of the Green 
function with respect to its two arguments |21j . The result is 

$(r 2 , a) = ^ / dSfo, B (r)V ■ ETAr, z = 0,t = s/c) 

2lT JSb 

/ dSVfaMr) ■E s ^{r,z = Q,t = s/c) 
2n Js n 

I dSV<l>2,B-V{<f> 1 ,A-<i>l,B)+ I dSVfo.B-V^B 



2: 

Qi 
~2tt 



Sb 
6( S ) 



(33) 



where in the second integral we integrated by parts and used the fact that 4>2.b vanishes on 
the wall of pipe B. Our final result for the wake becomes 



w = ^-6(s)I, 



where 



1=1 V<h,B{r)'V<h,B(r)dS- I V0 M (r)-V^ 2 , B (r)d5. 

JSb J Sap 

The impedance corresponding to this wake is 



(34) 



(35) 



(36) 



It is easy to see that this impedance is consistent with the symmetrized formula Eq. (251. 



VI. TRANSVERSE IMPEDANCE AND SMALL OFFSET OF PARTICLES 

Knowledge of the longitudinal impedance allows one to compute the transverse impedance 
using the Panofsky-Wenzel theorem [22 . In the general case, the transverse impedance is 
represented by a vector Z j_ perpendicular to the particle's orbit, and is given by 

Z_l = -V r ,Zu , (37) 

where V r2 is the operator "nabla" that differentiates with respect to the coordinates r 2 of 
the trailing particle, V r2 = xd/dx-i + yd/dy2, with x and y being the unit vectors in x 
and y directions 2 . 



In this paper we use the following definitions of the transverse wake w± 

c r°° 

wj_(ri. V2,s) = — / \Ei x( r 2, z = ct — s, t) + z X H±(r2, z = ct — s, t)] dt , 

91 J-oc 

and the transverse impedance Zj_(ri,r2,u>) = — (i/c) ds w± (ri, r2, s)e luJ3 / c . 



12 



In applications, it is typically assumed that there is a symmetry axis in the system and 
the beam has a small offset relative to this axis compared to the transverse size of the pipe. 
In this case, one can expand the impedance in Taylor series. The leading terms in the 
transverse impedance in this case are linear in offsets of the leading and trailing particles 3 . 
We will now derive expressions for the transverse impedance in this approximation. 

For small offsets of the leading and trailing particles we can expand the delta functions 
inEqs. (fl9) 



S(r - n.) w S(r) - r x ■ VS(r) + -(n • V) 2 c5(r) , 
5(r - r 2 ) sa S(r) - r 2 ■ VS(r) + ~(r 2 ■ V) 2 c5(r) , 



and represent each potential as a sum of a monopole part 
quadrupole part <// q ) , 



a dipole part 



0! = (t> {m) + 01 d) + 

where the parts satisfy the equations: 
V 2 (m) = -4nS(r) , 



j(<0 



(38) 



and a 



(39) 



/2^,(q) 



47rri • V<S(r) , 
-2tt(t-i • V) 2 <5(r) , 



V 



2 4 d) ^ 

72j,(q) 



47rr 2 • V(5(r) , 
= -27r(r 2 • V) 2 5(r) . 



(40) 



Representation (39 1 generates many terms in Eq. (36 1. The linear term in r 2 after sub- 



stitution into the Panofsky-Wenzel Eq. (37 1 gives rise to a transverse impedance which 



corresponds to a kick on particle 2 when both particles travel along the reference orbit with- 
out offset. We will call this impedance the transverse monopole impedance and denote the 
corresponding longitudinal impedance by Zu m , where 



Zn 



1 

2ttc 



v4 d i dS - 



1 

2vrc 



V0d 



V<4 d i dS . 



(41) 



In many practical cases the structure geometry has both up-down and right-left symmetry 
and the design orbit is on the symmetry line. In such a case the transverse monopole 
impedance vanishes. The longitudinal impedance which gives the usual transverse wake 
then has a term that is proportional to the product of vector components of r 2 and r\ and 
one that is quadratic in r 2 (no dependence on n). According to accepted terminology, we 
will call the former (in the case of impedance) the dipole component, Z^ d , and the latter 
the quadrupole component, Z\\^ q . The total impedance, the sum of the two components, is 
given by 



Z«=Z« d + Z n 



where 



Zn 



2ttc 



dS - 



2ttc 



V^ 2 %dS, 



(42) 



(43) 



As is well known, for axisymmetric systems the transverse impedance does not depend on the offset of 
the trailing particle. This however is not true for systems which are not axisymmetric. 
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and 



Sb 



'l.B 



V^ B dS 



2nc 



(44) 



Note that for the special case of cylindrically symmetric structures the quadrupole compo- 
nent is identically equal to zero. 



Equations (|3T[) , (41 1, (43 1, and (44) allows one to calculate all the component of the 



transverse impedance for a transition of arbitrary geometry. 



VII. EXAMPLE IMPEDANCE CALCULATIONS IN THE OPTICAL REGIME 



In this section we will show how to calculate the longitudinal and transverse impedances 
for several simple cases of axisymmetric systems using results of Section |VI| 

We first calculate the longitudinal impedance for a short, round collimator shown in 
Fig. with pipe radius a and collimator radius b (b < a). In this case both charges are 




a) 



b) 




FIG. 4: A short collimator (a), a step-in (b), and a step-out (c) transition in axisymmetric geometry. 
The particles move in the +z direction. For the collimator a denotes the pipe radius and b the 
collimator radius, with b < a. For the step transitions, a denotes the radius of the incoming pipe 
and b the radius of the outgoing pipe, with b < a for the step-in and b > a for the step-out transition. 



located on the axis of the pipe, and the solution to Eqs. (19), in cylindrical coordinates, is 

r 

(45) 



b = -21n- 
a 



Substituting this solution into Eq. ( 36 ) yields 



z ^L! s ™ 2ds ~LS s ™ 2dS 

= ^-l (V^) 2 ^=*ln£. 

2™Js B -s* P c b 



(46) 



This is a well known result for the collimator impedance in the high frequency limit. 
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The longitudinal impedance for step-in and step-out transitions shown in Fig. [4J) and 
|4p can easily be obtained from the previous equations. In both cases the potential <j> is 
— 21n(r/a) and — 21n(r/6) for pipes of radius a and b, respectively. However, for the step-in 
transition, the cross-section Sb coincides with S ap , and the two integrals in Eq. (35 1 cancel, 



giving a total impedance of zero. For a step-out transition, the difference of these integrals 
is the same as given by the first line in Eq. (46 1, and the impedance is equal to that of a 



collimator with the aperture equal to the radius of pipe B. Again, both these results are 
known in the literature, and we have derived them here to demonstrate that the optical 
approximation agrees with previously obtained results. 

For the transverse impedance of a round collimator and a step-in and a step-out transition 
we first note that only the dipole term in Eq. (42 1 contributes to the impedance. The 



quadrupole term in this equation vanishes because the monopole potential 4>( m ^ in Eq. ( 44 1 
has no angular dependence while the quadrupole potential </>( q ' has an angular dependence 
oc cos 29 (here 9 is the azimuthal angle in cylindrical coordinates) . 

For a round collimator, pipe A and pipe B have the same radius a, hence 



9 l.A 



(d) 
l.B 



and 



This reduces the integration in Eq. ( 43 1 to one over the collimator area 



1 

2ttc 



SB — S a -r 



(47) 



We assume now that both leading and trailing particles are offset in the direction of the x 
axis and define function tp such that 



The function ip A satisfies the following equation 

V 2 Vm = AnS'(x)8(y) , 



(48) 



(49) 



where the prime denotes derivative with respect to the argument, with the boundary con- 
dition ipA — at r = a. The solution is 



ipA = —2x 



1 



1 



x 2 + y 2 



= — 2rcos( 



1 



1 



(50) 



The expressions ( 48 1 are next substituted into Eq. ( 47 1 and integrated over the cross 



sectional area of the collimator. The calculation can be simplified if one uses the identity 



VipA-y^AdS = V ■ (i/) A VipA)dS - tpA^ ^AdS. (51) 

SB~S a p JSb — S ap J Sb — S ap 

The second integral on the right hand side vanishes because V 2 V>a = in any region that 
does not include the point x = y = 0, and the first one can be cast into an integral over the 
edge of the collimator r = b 



V-{^ A V^ A )dS = -b [ \ A ^±d6=^r (l h 



'Sb — Sap •* 

which gives for the longitudinal dipole impedance 



dr 



b 2 



7 2 ^ 



(52) 



(53) 
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Using now the Panofsky-Wenzel relation Eq. ( 37 1 we find the x component of the transverse 
impedance 

Z_L. X _ 2 

Zap 



Xi 



(54) 



The right hand side of this equations defines the transverse impedance per unit offset X\- 



this quantity is traditionally referred to as the transverse impedance. The result Eq. (54 1 
agrees with Ref. [12] . 

For the step-in and step-out transitions, Eqs. (48 1 and (50 1 define the dipole potentials 
for pipe A; the potentials for pipe B are obtained from these formulae by exchanging b and 
a, 



j(d) 

l,B 



Xi_1p 



B ■ 



6 (d) 
°2.B 



-2x 



IP 



X24>B 

= — 2r cos ( 



(55) 



Note that the difference Aip = ipA 
direction 



ipB corresponds to a uniform electric field in the x 



Aip = 2x 



1 

P 



(56) 



For the step-in transition S ap 
reduced to 



Zn 



Sb, and after a simple transformation Eq. (43 1 can be 
f (VAV>- Vtp B ) dS. 

JSn 



X\Xi 
2lTC 



(57) 



Using the same integration by parts as in Eq. (51 1 we obtain 

I (VAt/> ■ Vtp B ) dS = I V ■ (AipVip B ) dS- I (Ai/)V 2 xp B ) dS 

J Sb ** S j3 J Sb 

dips 



= b 



Aip - 



Or 



r—b 



de-Aii / Aip5' {x)8{y)dS 



Sb 



= 0, (58) 

which means that the longitudinal dipole impedance, and hence the transverse impedance, 
is equal to zero for the step-in transition. We remind the reader that the longitudinal 
impedance for this case also vanishes. 

For the step-out transition 5 ap = Sa- We represent the integral over S ap in Eq. (43 1 
as the difference of the integrals over Sb and the integral over Sb — S ap . Then the first 
term in this equation combined with the integrals over Sb gives zero, as follows from the 
calculations for the step-in transition. Hence 



Z \\.d = 



X\X 2 

2wc 

X\X 2 

2ttc 

X\X 2 ^ 
2lTC 

4:XlX 2 



(VtpA ■ V^ B ) dS 



Sb — Sz 



Sb — Ss. 
2tt 



V • (iPaV^b) dS 



TP. 



dipt 



Or 



(W 



r—b 



l 

Jp 



(59) 
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For the transverse impedance in this case we find 



4 / a 



which agrees with the result of Ref . |5 . 



b 2 



(60) 



VIII. PILLBOX CAVITY AND RELATION BETWEEN OPTICAL AND 
DIFFRACTION REGIMES 



Let us consider now the pillbox cavity shown in Fig. [5] Diffraction theory gives the 



A 



It 



FIG. 5: An axisymmetric pillbox cavity, 
longitudinal impedance for a cavity as (see, e.g. |19j ) 



, diffraction 



2(1 + $) / I 



,1/2 



cb 2 uj 



(61) 



where b is the pipe radius and / is the length of the cavity. 

The optical theory predicts zero impedance for the pillbox cavity. Indeed, in this case, 
all the three cross-sections Sa, S ap and Sb are equal (regarding the determination of S ap 
in this case see the footnote on page|9|, and Eq. (35 1 immediately gives a zero result. The 
reason for the optical approximation not reproducing the result of the diffraction theory is 
that Eq. (61 1 corresponds to the next order approximation in the small parameter Al/b 2 , 



which is beyond the applicability limit of the optical regime. Indeed, if we take the ratio of 
the impedance Eq. (61 1 to a typical optical impedance Z|| i0ptica i ~ l/ c [ see , s.g. Eq. (46 1], 
we obtain 



(62) 




Hence the parameter y/lX/b 2 indicates the accuracy of the optical approximation: one can 
expect that corrections to the optical regime due to diffraction effects will be on the order 
of this parameter 4 . 

Note that one can find in the literature references to the impedance of the collimator 
given by Eq. ( 46 1 as a diffraction impedance (the authors of this paper have also used 
this terminology in the past). The terminology introduced in this paper distinguishes the 
optical regime from the diffraction regime, and we believe that it better describes the physics 
involved. 



4 For the geometries considered in Section |VII| the length parameter £ should be denned as I ~ b, and the 
accuracy of the optical approximation is on the order of y/X/b. 
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IX. CONCLUSION 



In this paper we have introduced an optical approximation into the theory of impedance 
calculation, valid in the limit of high frequencies. This approximation neglects diffraction 
effects in the radiation process, and is conceptually equivalent to the approximation of 
geometric optics in electromagnetic theory. Using this approximation, we have derived 
equations for the longitudinal impedance for arbitrary offsets of the source and test particles 
with respect to a reference orbit. With the help of the Panofsky-Wenzel theorem we have 
also obtained expressions for the transverse impedance (also for arbitrary offsets). We 
further simplified these expressions for the case of the small offsets that are typical for 
practical applications. Our final expressions for the impedance, in the general case, involve 
two dimensional integrations over various cross-sections of the transition. 

We have, in addition, demonstrated for several simple examples how our method is applied 
to the calculation of impedances for simple axisymmetric geometries that have been studied 
in the past. Finally, we discussed the accuracy of the optical approximation and its relation 
to the diffraction regime in the theory of impedance. 
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